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A wide range of thermoreversible gels are prepared by cooling down to ambient temperature hot aqueous 
polymer solutions. During the sol-gel transition, such materials may experience a volume contraction which is 
traditionally overlooked as rheological measurements are usually performed in geometries of constant volume. 
In this article, we revisit the formation of 1.5% wt. agar gels through a series of benchmark rheological 
experiments performed with a plate-plate geometry. We demonstrate on that particular gel of polysaccharides 
that the contraction associated withe the sol/gel transition cannot be neglected. Indeed, imposing a constant 
gap width during the gelation results in the strain hardening of the sample, as evidenced by the large negative 
normal force that develops. Such hardening leads to the slow drift in time of the gel elastic modulus G 
towards ever larger values, and thus to an erroneous estimate of G. As an alternative, we show that imposing 
a constant normal force equals to zero during the gelation, instead of a constant gap width, suppresses the 
hardening as the decrease of the gap compensates for the sample contraction. As such, imposing a zero normal 
force is a more reliable way to measure the linear properties of agar gels, which we prove to work equally well 
with rough and smooth boundary conditions. Using normal force controlled rheology, we then investigate the 
impact of thermal history on 1.5% wt. agar gels. We show that neither the value of the cooling rate, nor the 
introduction of a constant temperature stage during the cooling process influence the gel elastic properties. 
Instead, G only depends on the terminal temperature reached at the end of the cooling ramp, as confirmed 
by direct imaging of the gel microstructure by cryoelectron microscopy. Finally, we also discuss two subtle 
artifacts associated with the use of duralumin plates that may interfere with the rheological measurements 
of agar gelation. We show that (*) the corrosion of duralumin by the aqueous solution, and ( ii ) the slow 
migration of the oil rim added around the sample to prevent evaporation, may both lead separately to a 
premature and artificial growth of G that should not be misinterpreted as the formation of a pre-gel. The 
present work offers an extensive review of the technical difficulties associated with the rheology of hydrogels 
and paves the way for a systematic use of normal force controlled rheology to monitor non-isochoric processes. 


I. INTRODUCTION 

Bulk rheology is a powerful tool to monitor the linear 
properties of complex fluids such as polymer and sur¬ 
factant solutions, emulsions, colloidal suspensions, gels, 
etc!- ; . Although most traditional rheology experiments 
are performed at constant temperature, several processes 
involve large temperature variations that make rheologi¬ 
cal measurements much more challenging. A textbook 
case of rheological experiments involving temperature 
variations, is the sol-gel transition. Popular examples 
include the jamming transition in suspensions of ther- 
m osen sitive particleiP"^], the precipitation of waxy crude 
oil® 7 and the heat-induced formation of numerous gels 
such as biop olyme r gels^OS gels made of oil-in-water 
nanoemulsionskLU^ e tc. In all of these examples, gela¬ 
tion is traditionally monitored in a geometry of constant 
gap with roughened boundary conditions to prevent wall 
slip. Small amplitude oscillatory shear (SAOS) is used 
to measure the linear elastic and viscous moduli, respec¬ 
tively G and G ", and the thermal dilation of the shear 
cell due to the temperature variation is taken into ac- 
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count through a gap compensation procedure that is now 
commonly implemented on rheometers. 

Nevertheless, a key technical assumption in such mea¬ 
surements is that the contraction of the sample during 
the sol-gel (or the liquid-solid) transition is negligible. 
Such assumption is not always valid as phase transition 
often comes with a volume contraction of the sample. In 
the latter case, if the phase transition is monitored in a 
geometry of constant volume, the sample will experience 
strain hardening as it deforms and remains in contact 
with the cell boundaries. In a more subtle way, the con¬ 
traction may also lead either to the complete debonding 
of the sample from the walls of the shear cell as observed 
for very rigid samples, or to the formation of lubrication 
patches at the gel/geometry interface for soft hydrogels 13 
which in both cases entail erroneous measurements of 
the sample elastic properties. Such issue is particularly 
important for biopolymer gels made of polysacharides 
and/or proteins which are the topic of the present work 
as they experience volume variations during the sol-gel 
transition, and are prone to solvent release under exter¬ 
nal stress. 

Quantitative measurements of the linear propert ies of 
such soft solids are crucial in numerous application d 14 l 15 1 
For instance, biopolymer gels serve as substrate for bac¬ 
terial growth and cell culture. The gel stiffness strongly 
impacts the cell locomotion, a phenomenon that is called 
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“durotaxis”^, as well as stems cell differentiation 17 
which requires reliable methods to quantify the elastic 
properties of gels, especially when the latter is prepared 
through a complex thermal history. To our knowledge, 
most if not all the gelation experiments on thermore- 
versible gels reported in the literature have been per¬ 
formed at constant gap width, either in parallel-plate 
and_cone-and-plate geometrj®HISl or j n a Taylor-Couette 
cell™] None of these experiments takes into account 
the contraction of the sample during gelation, although 
the loss of adhesion between the gel and the wall of the 
shear cell has been already reported as a major issue 1 ^. 
The use of sand paper, although widespread, do not com¬ 
pensate for the sample contraction. Moreover, the pres¬ 
ence of such a medium sandwiched between the Peltier 
plate and the sample further impairs the temperature 
control of the sample, while the glue attaching the sand¬ 
paper poorly withstands high temperatures. As a result 
of the sample contraction, measurements of the viscoelas¬ 
tic moduli in different geometries are often inconsistent 
which urges experimentalists to use alternative methods 
to monitor more accurately non-isochoric processes such 
as the gelation of biopolymer gels. 

The purpose of the present paper is twofold. First, it 
aims to illustrate in a series of benchmark experiments in 
the same spirit as a recent work by Ewoldt et alP3 that 
the use of a normal force controlled procedure is more 
adapted to monitor the gelation of thermoreversible gels 
than the traditional constant gap protocol. Normal force 
controlled rheology has been applied recently with suc¬ 
cess at constant temperature to the study of dense gran¬ 
ular suspensions^SHMl. However, if such method is also 
advertised by Anton Paar to probe soft solids in “World 
of Rheology”^, or appears in the TA instruments help 
files, the present manuscript is to our knowledge among 
the first quantitative report regarding normal force con¬ 
trolled rheology applied to thermoreversible gels. Here 
we show on a 1.5% wt. agar gel, that applying a constant 
normal force equals to zero to the sample enclosed in a 
parallel-plate geometry prevents any strain hardening or 
debonding of the gel from the plates, thus providing pre¬ 
cise measurements of both the gelation dynamics and the 
gel elastic properties. Furthermore, we show that such 
method works equally well with both rough and smooth 
boundary conditions, and that as such it can be used 
advantageously with bare metallic surfaces without any 
sandpaper. Second, we build upon the first technical part 
and use the zero normal force protocol to study the im¬ 
pact of the thermal history of agar sols on the mechanical 
properties of the subsequent gels. We show that the cool¬ 
ing rate of the solution sets the gelation temperature and 
influences the gel contraction, without affecting signifi¬ 
cantly the final microstructure of the gel. Furthermore, 
the steady-state value of the gel elastic modulus is shown 
to be determined only by the final temperature reached 
at the end of the cooling ramp, and to be insensitive to 
the introduction of a constant temperature stage during 
the colling process. 


Finally, we highlight two artifacts that may take place 
during the gelation dynamics of agar gels, whether the 
gelation is monitored at constant gap or at constant nor¬ 
mal force. First, we show that the slow oxidation of the 
metallic walls of the shear cell by the agar sol may lead to 
an artificial increase of G' while the sample is still liquid. 
Second, we demonstrate that the oil rim surrounding the 
sample that is traditionally used to prevent the solvent 
evaporation, may slowly invade the gap also triggering 
a premature increase of the gel elastic modulus before 
the actual start of the gelation. Such an artifact leads 
to an erroneous estimate of the crossing of G' and G", 
but hardly impacts the measurement of the gel final elas¬ 
tic properties, as long as the terminal elastic modulus is 
larger than a few hundred pascals. 

The paper is organized as follows. Section [TT| describes 
the agar samples, the experimental setup and the zero 
normal force procedure. We also discuss the effect of 
the strain amplitude and the impact of boundary con¬ 
ditions on the rheological measurements. In section m 
the zero normal force protocol is applied to determine the 
impact of thermal history on the mechanical properties 
of agar gels. Finally, section |TV| provides a short discus¬ 
sion of these results and emphasizes future applications 
of normal-force controlled rheology. 


II. MATERIAL AND METHODS 
A. Agar gels 

Agar is extracted from a marine red algae and consists 
in a mixture of polysaccharides among which the gelling 
component is agarose, a polymer composed of two dif¬ 
ferent galactose subunitj 32 ! 33 ! When agar is dissolved in 
boiling water, the agarose is responsible for the forma¬ 
tion of a crosslinked network upon cooling. The gelation 
involves a liquid-liquid phase separation known as spin- 
odal decomposition that corresponds to the formation of 
polymer-rich and solvent-rich regionsj^^H, together with 
a conformational change of the a garose molecules which 
self-associate via hydrogen bonds 37 2H. The exact for¬ 
mation scenario of the gel results from a subtle compe¬ 
tition between these three distinct processed^ which in 
turn strongly depends on the agarose con centrati on, its 
molecular weight and the thermal histord ^ 5 ! 40 ! 41 1. The 
spinodal demixion is favored for agarose concentrations 
lower than 2% wt., whereas the gelation occurs through 
a more direct scenario at larger concentrations. Further¬ 
more, although the exact conformation of the agarose 
molecules inside the gel network is still debated (single 
vs. double helices}® S3, the microstructure of agar gels 
at a coarser scale is consensually described as a p orous 
network made of bundles of agarose chain4® 48 * 49 l. The 
latter microstructure is filled with water and displays 
a brittle-like mechanical behavior including the forma¬ 
tion of macroscopic fractures at large enough straind^HUl. 
Agar gels are also prone to release water at restful, or 
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FIG. 1. (Color online) Temporal evolution of the normal force Fjy and the gap width e (a)-(c) together with the elastic and 
viscous moduli (d)-(e) during the gelation of a 1.5% wt. agar solution induced by decreasing the temperature from T = 70° C 
to 20°C at a cooling rate T = l°C/min. The gelation experiment is repeated three times with different protocols. The 
first column illustrates the gelation performed without any particular precaution. Although the gap width is set initially at 
eo = 500 (im at T = 20°C, the true gap width varies by 65 pm due to the thermal dilation of the geometry, resulting in an 
lower terminal value of the elastic modulus G'f = 11.7 kPa. The second column shows a gelation experiment during which the 
temperature compensation mode is active and compensates for the thermal dilation of the geometry alone. The gap width is 
truly constant during the whole gelation (eo = 500 pm). Yet, the normal force is still negative which proves that the sample 
contracts during the gelation. Such contraction leads to strain hardening as evidenced by the continuous increase of the elastic 
modulus ( G'f ~ 28.0 kPa at t = 5 h). Finally, the third column illustrates a gelation experiment during which a constant 
normal force is applied to the sample Fn = (0.0 ± 0.1) N, while the temperature compensation mode is active. The gap width 
of initial value eo = 500 pm, decreases by 1% which compensates for the sample contraction (G'f = 25.9 kPa). For the three 
experiments, the strain amplitude is 7 = 0.01 %, the frequency is / = 1 Hz and the solvent trap is filled with water. 


under external deformation as low as a few pourcents— 1 
which makes their rheological study particularly chal¬ 
lenging and motivates the present study. 

Here we focus on hydrogels made with 1.5 % wt of 
agar composed at 70% of agarose (BioMerieux). Agar 
is added as a powder to warm deionized water, and the 
mixture is brought to a boil whilst stirring mildly for 
about 10 minutes. The solution is then cooled down to 
T = 80° C and kept at this temperature up to three days 
during which several samples are taken out to prepare 
gels inside the gap of a rheometer, as detailed below. 
After three days at 80° C, the polymer solution starts 
suffering chemical aging as e videnc ed by the formation 
of gels of lower elastic modulf^I 57 . Fresh solutions are 
therefore prepared regularly, and all the data reported 
here have been obtained from fresh batches, younger than 
three days. 


B. Rheological setup 

Rheological measurements are conducted in a plate- 
plate geometry of diameter 40 mm driven by a stress- 
controlled rheometer (DHR-2, TA Instruments). The 


upper moving plate made of duralumin is sand-blasted 
and displays a surface roughness of 4 ± 2 /xm, measured 
by Atomic Force Microscopy. The bottom plate consists 
in a smooth and Teflon-coated Peltier unit that allows 
us to control the temperature of the sample. Evapo¬ 
ration is minimized either by using a solvent trap con¬ 
taining deionized water and placed on top of the upper 
plate, or by adding a thin layer of sunflower seed oil from 
Helianthus annuus (Sigma Aldrich) around the sample 
depending on the duration of the experiment (see sec¬ 
tion IIF for more details). Gelation experiments are con¬ 


ducted as follows: the agar solution is first introduced at 
T = 80° C in the gap of the plate-plate geometry that 
has been preheated at T = 70°C. The temperature is 
then decreased at a constant rate T = l°C/min down 
to Tf = 20° C (unless stated otherwise) and maintained 
at this temperature up to several hours to make sure 
the gelation is complete. In the meantime, oscillations of 
small amplitude performed at a frequency / = 1 Hz allow 
us to monitor the temporal evolution of the gel viscoelas¬ 
tic properties. Gelation experiments are conducted either 
while imposing a constant gap width (e = 500 /xrn), or a 
constant normal force (Fjv = 0.0 ± 0.1 N) while the gap 
width is thus free to vary with the gel thickness. Finally, 
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a few tests are performed on deionized water at constant 
temperature using either one or two transparent plates 
made of PMMA. In that case, images of the sample are 
taken from below with a webcam (Logitech HD c920) by 
means of a flat mirror placed at an angle of 45° with 
respect to the horizontal plates. 


C. Gap compensation and zero normal force protocol 

Let us start by illustrating the impact of the thermal 
expansion of the shear cell on the rheological measure¬ 
ments. For this first experiment, the gap width is set 
to e = 500 /tm at T = 20°C. However, due to the ther¬ 
mal dilation of the plates, the gap width decreases by 
1.3 /jm/°C for increasing temperature, as determined 
by calibration. The gap width is therefore 435 /im at 
the moment the polymer solution is introduced inside 
the pre-heated geometry. Once the polymer solution is 
loaded, the temperature is decreased at a constant rate 
T = l°C/min down to Tf = 20°C while we monitor the 
evolution of the normal force [Fig. [l|a)] and apply oscil¬ 
lations of small amplitude with a strain 7 = 0.01 % to de¬ 
termine the elastic and viscous moduli of the gel, respec¬ 
tively G' and G" [Fig. Qd)]. During the first half hour, 
while T > 33° C the agar solution remains liquid, and 
the normal force F N is close to zero as the gap increases 
back to the initial value of 500 /im. At i ~ 0.6 h, the gel 
starts forming as evidenced by the growth of the elastic 
modulus which becomes larger than the viscous modulus. 
Concomitantly, the normal force Fjv shows negative val¬ 
ues indicating that the gel is pulling down on the upper 
plate, until a strong upturn occurs at t = 0.75 h before 
the end of the cooling phase (t = 0.8 h), indicating a bru¬ 
tal change in the contact between the gel and the plates, 
as discussed below. Moreover, G' goes through a max¬ 
imum at t = 0.88 h after the temperature has reached 
Tf = 20° C, and relaxes towards a steady state value 
G'f = 11.7 kPa, while Fn remains negative. The forma¬ 
tion of the gel appears as over after 2.5 h, and the gel 
displays a solid-like behavior as evidenced by G' G". 

In fact, the above procedure which does not compen¬ 
sate for the thermal dilation of the geometry leads to 
artificially low estimates of G'. Indeed, we have repeated 
the gelation experiment taking into account the variation 
of the gap thickness due to thermal dilation [Fig. [ljb)] . 
In that case, later referred to as the “temperature com¬ 
pensation mode”, the gap width remains constant during 
the entire experiment and the terminal value of the elastic 
modulus is G'f = 28.0 kPa [Fig.jlje)], instead of 11.7 kPa 
in the absence of compensation ior the thermal dilation 
of the geometry [Fig. |T](d)] . This experiment illustrates 
how dramatic and misleading can be the effect of the 
thermal dilation of the geometry on the measurements of 
the gel linear properties. As a consequence, the abrupt 
upturn of the normal force observed in Fig. Qa) in the 
absence of any compensation results from a partial loss 
of contact between the sample and the geometry, and/or 


the formation of lubrication spots between the gel and 
the plate due to the large values of F,v |. Indeed, from 
the maximum value of |Fjv| ~ 8 N in Fig. [lja) , one can 
estimate the deformation experienced by the gel with an 
elastic modulus of 28 kPa to be of about 20% which is 
large enough to trigger the release of solvent and the par¬ 
tial lubrication of the gel-plate interface. 

Nonetheless, the temperature compensation mode does 
not lead to fully satisfying rheological measurements. In¬ 
deed, one observes that the normal force still exhibits 
negative values during the gelation although we compen¬ 
sate for the thermal dilation of the geometry [Fig. [ljb)]. 
This result shows that the gel contracts during the sol- 
gel transition and pulls down on the upper plate. As a 
result, G' displays a weak logarithmic increase even after 
F/v has reached a plateau [Fig. [lje)] which is the signa¬ 
ture of the strain hardening of the agar gel induced by 
the sample contraction [see Fig.[l3]in the appendix for a 
supplemental experiment and the corresponding discus¬ 
sion] . Such a slow drift of G' towards ever larger values is 
also visible in other studies conducted on gels containing 
agar(ose[P^®], but usually unnoticed and not discussed 
as rheological data are most often plotted in semiloga- 
rithmic scale which flattens G' artificially. 

To compensate the effect of both the thermal expan¬ 
sion of the geometry and the contraction of the sam¬ 
ple during the phase transition, the gel formation is 
now monitored with the temperature compensation mode 
switched on, while imposing a controlled normal force 
F/v = (0.0 ±0.1) N instead of a constant gap width. The 
gap width, which initial value is set to eo = 500 /tm may 
now vary during the experiment and indeed, we observe 
a gap decrease of 4 /im concomitantly to the growth of 
G' [Fig. [ljc)]. Such a gap decrease of about 1% allows 
the upper plate to stay in contact with the gel while the 
latter contracts during the gelation. Furthermore, the 
elastic modulus does not drift anymore and reaches a 
constant value G/ = 25.9 kPa at the end of the cooling 
phase [Fig. [lj'f)] , that is smaller than the one recorded 
in the constant gap experiment [Fig. |TJ[e)] . The latter 
observation also confirms the hardening scenario invoked 
in Fig. |TJ(e) . The error bar associated with G'f is mainly 
given by the ability of the experimentalist to load the 
shear cell with the same amount of hot agar solution. 
Repeating the same experiment about ten times on sam¬ 
ples extracted from different batches leads to compatible 
results within 10%. There is thus only a slight difference 
between the value of G'f measured with a constant gap 
(after 5 h), and the value measured with a constant nor¬ 
mal forced However the value of G/ measured with a 
constant gap width protocol is a function of time and as 
such depends on the patience of the experimentalist. In 
conclusion, monitoring the gelation while imposing a con¬ 
stant normal force equals to zero prevents the strain hard¬ 
ening of the gel and appears as a more accurate method 
to determine the gel elastic properties. This method, that 
is further referred to as “the zero normal force protocol” 
is the one that is used in the rest of the manuscript. 
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FIG. 2. (Color online) The first column shows the temporal evolution of (a) the gap width e while imposing a constant normal 
Fn = (0.0 ± 0.1) N during the gelation of a 1.5% wt. agar solution induced by decreasing the temperature from T = 70°C to 
20°C at T = l°C/min (relative decrease of the gap width Ae/e = 1.4%). (b) Evolution of the elastic and viscous moduli in 
semilogarithmic scale determined through small amplitude oscillations with 7 = 0.01 % and f = 1 Hz (Crossover temperature 
of G' and G": T g = 35.0°C; terminal value of the elastic modulus: G'j = 21.0 kPa). The second column shows the same 
experiment performed at / = 1 Hz, but during which the strain 7 is adapted to the value of the elastic modulus for a better 
resolution from the early stage of the experiment while the sample is still liquid: the strain amplitude is 7 = 1 % for G' < 1 Pa, 
7 = 0.1% for 1 Pa < G' < 10 Pa, and 7 = 0.01 % for G' > 10 Pa (Ae/e = 0.6 %; T g = 35.2°C; G/ = 21.6 kPa). For both 
experiments, the solvent trap is filled with water. 


In the rest of section[TT]we show that the terminal value 
of G' determined with the zero normal force protocol 
is robust and independent of the surface roughness of 
the plates. We also discuss in depth several artifacts 
related to the strain amplitude applied to measure the 
viscoelastic moduli, and to physico-chemical phenomena 


related to the boundary conditions. Finally, in section III 


we use the zero normal force protocol to investigate the 
impact of thermal history on the linear properties of a 
1.5% wt. agar gel. 


D. Strain adapted protocol 


Here we briefly discuss the impact of the strain ampli¬ 
tude 7 applied to monitor the gel formation. Fig. [2^b) 
displays in a semi-logarithmic plot the temporal evolu¬ 
tion of the elastic and viscous moduli during a gelation 
conducted with the zero normal force protocol and deter¬ 
mined with 7 = 0.01 % and / = 1 Hz [same protocol as in 
Fig.Qf)]. Although the sample is liquid-like during the 
first half hour, such a small deformation leads to noisy 
and comparable values of G' and G" until both mod¬ 
uli exceed about 1 Pa. This experimental obse rvat ion, 
which is commonly encountered in the literature 63 ^, is 
a consequence of the too small applied value of the strain 
amplitude. Indeed, the minimum measurable value of G' 


is inversely proportional to 7 and given by the following 
expression: G% n = 4 T min /( 37 rI? 3 7 ), where R denotes the 
plate radius and T m - ln the smallest measurable torque 2 '. 
Here, T m i n = 2 nN.m (for a DHR-2) and with 7 = 0.01 % 
we find G% n = 1 Pa which is in excellent agreement with 
the minimal significant stress value observed in Fig.[2](b). 
Imposing a larger strain before the start of the gelation, 
such as 7 = 1 % allows us to decrease G' min down to 
0.01 Pa and measure consistent values for the viscoelas¬ 
tic moduli during the first half hour, i.e. G' = 0 since the 
sample is liquid [Fig. [2jd)]. The strain amplitude is then 
decreased when G' increases as the gelation starts, so as 
to remain in the linear regime and prevent any strain- 
induced release of solvent and/or debonding of the sam¬ 
ple from the geometry. Indeed, strain sweep experiments 
performed on gels which have been prepared under os¬ 
cillations of different amplitudes 7 indicate that strain 
amplitudes as low as 7 = 0 . 01 % are necessary for an opti¬ 


mal adhesion between the gel and the plates (see Fig. 14 


and related discussion in the appendix which proves that 
the gel adhesion to the plates decreases for increasing 
strain values 7 applied during the gelation). As a conse¬ 
quence, the following protocol leads to reproducible mea¬ 
surements: 7 = 1 % for G' < 1 Pa, 7 = 0.1 % for 
1 Pa< G' < 10 Pa and 7 = 0.01 % for G' > 10 Pa, and 
confirms that as long as the solution is liquid, only G" 
displays non-negligible values. After gelation, the termi- 
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FIG. 3. (Color online) (a) Picture of the duralumin rotor which has been polished with a suspension of ultrafine aluminium 
oxide particles. The scale is fixed by the plate diameter which is 40 mm. (b) Evolution of the viscoelastic moduli G' and G" 
vs. time during the gelation of a 1.5% wt. agar solution induced by decreasing the temperature from T = 70°C to 20°C at 
T = l°C/min, under controlled normal force Fjv = (0.0 ± 0.1) N (initial gap size eo = 500 pm). The experiment is performed 
twice, each time with symmetric boundary conditions: once with smooth walls, i.e. a polished duralumin upper plate pictured 
in (a) and a Teflon coated bottom plate (red and blue symbols) and once with rough walls, i.e. sand-blasted duralumin plates 
(black and gray symbols). In both cases the relative gap decrease is 1%, the crossover temperature of G' and G" is T g = 35.7°C, 
and the terminal values of the elastic modulus are respectively G'f = 23.8 kPa and 23.6 kPa. Both experiments are perfomed 
with a solvent trap filled with water. 


nal value of G' obtained with the strain adapted protocol 
is compatible within error bars with the value measured 
with a constant strain 7 = 0.01 % [Fig. [ljf )]. Therefore 
adapting the value of strain amplitude during the gela¬ 
tion offers the opportunity to determine more accurately 
the crossover between G' and G" which is often used as 
a good estimate of the gelation poinlP^ and which occurs 
here at T g = 35.2° C. Such a strain-adapted protocol is 
systematically applied during the experiments reported 
in the rest of the manuscript. 

E. Boundary conditions and related artifacts 
1. Surface roughness 

Let us now discuss the impact of the surface roughness 
of the plates. The experiments reported in Fig. [I] and [2] 
have been performed with mixed boundary conditions 
(BC), namely a rough metallic upper plate and a smooth 
(Teflon-coated) bottom Peltier plate. To further test the 
robustness of the zero normal force protocol, we have per¬ 
formed gelation experiments with symmetric boundary 
conditions. On the one hand, experiments with symmet¬ 
rically smooth boundary conditions are performed by re¬ 
placing the sand-blasted duralumin rotor with a duralu¬ 
min plate polished on a rotating wheel with a suspension 
of ultrafine aluminum oxide particles to produce a mir¬ 
ror like surface pictured in Fig.[3|a) and of typical rough¬ 
ness 10±15 nm, as determined from AFM measurements. 
The bottom plate is the smooth and Teflon-coated Peltier 
plate used by default. On the other hand, experiments 


with symmetrically rough boundary conditions are per¬ 
formed by topping the Peltier plate with a sandblasted 
duralumin cover that ensures a good thermal conductiv¬ 
ity between the Peltier plate and the sample, and allows 
us to achieve a roughness of about 4 /im for the bottom 
plate similar to that of the upper sandblasted duralu¬ 
min plate. Both experiments reported in Fig. [3](b) lead 
to compatible values of G^ within error bars (rough BC: 
Gj = 23.6 kPa; smooth BC: G'f = 23.8 kPa) that are also 
in quantitative agreement with the elastic modulus deter¬ 
mined with mixed boundary conditions [see Fig. 0 f)]- 
These experiments demonstrate that the zero normal 
force protocol together with oscillations of small ampli¬ 
tude adapted to the value of the elastic modulus allows 
one to monitor the gelation dynamics in a plate-plate 
geometry with either smooth or sand-blasted boundary 
conditions. The zero normal force protocol thus appears 
as much simpler than the traditional constant gap width 
protocol which requires the use of sandpaper. In par¬ 
ticular, sandpaper does not compensate the sample con¬ 
traction during gelation and is a poor thermal conduc¬ 
tor. Moreover, it often displays adhesion issues to the 
Peltier plate at high temperature, hence inducing large 
error bars. 


2. Chemical properties of the plates surfaces 

This subsection concerns a subtle artifact associated 
with the use of duralumin plates. We have observed that 
gelations performed with a duralumin rotor which has 
been freshly sandblasted before the experiment, system- 
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FIG. 4. (Color online) (a) Temporal evolution of G' and G" during the gelation of a 1.5% wt. agar solution induced 
by decreasing the temperature from T — 70°C to 20°C at a cooling rate T = l°C/min and under controlled normal force 
Fn = (0.0 ±0.1) N (initial gap size eo = 500 /.tm, relative decrease of the gap width: Ae/e = 1%); terminal value of the elastic 
modulus: G'j = 23.1 kPa. The duralumin rotor has been sandblasted just before the experiment, and the gelation is performed 
with a solvent trap filled with water, (b) Picture of the sandblasted rotor after a few successive gelation experiments. Notice 
the iridescent regions at the periphery of the plate which corresponds to the oxidized area. The scale is set by the rotor diameter 
of 40 mm. 


atically show a premature increase of the elastic modulus. 
An example is pictured in Fig. Qa): G' increases from 
t ~ 0.1 h instead of t = 0.55 h [compare with Fig. [2][d)] , 
and tends toward about 2 Pa before following the steep 
increase previously associated with the gelation. Such a 
premature increase of G' progressively disappears when 
the experiment is repeated four or five times with the very 
same geometry, while the upper plate takes an iridescent 
appearance [Fig. |4](b)] . If the rotor is sandblasted anew, 
the premature increase of G' is again visible in the sub¬ 
sequent first couple of experiments. We attribute such a 
premature increase of G' to the oxidation of the upper 
duralumin plate by the agar aqueous solution. Nanopar¬ 
ticles of aluminium oxide detach from the surface of the 
upper plate to be advected through the sample toward 
the air/liquid interface, which increases the elasticity of 
the contact line and leads to the premature increase of 
G'. In this framework, the successive gelation experi¬ 
ments result in the passivation of the upper plate which 
strongly slows down any further oxidation of the geome¬ 
try and explains why the premature increase of G' is no 
longer observed in the subsequent experiments. 

To support this interpretation, we have performed the 
following test: the rotor is sandblasted anew, and the 
preheated gap is filled with distilled water instead of the 
agar sol. The temperature is decreased as usual from 
T = 70°C down to 20°C at a cooling rate T = l°C/min 
under a constant normal force Fn =(0.0±0.1) N and os¬ 
cillations of small amplitude (7 = 1 % and f = 1 Hz). Al¬ 
though no gelation can take place, one can see in Fig.[5](a) 
that after 5 min, the elastic modulus increases up to a few 
Pa, and remains larger than G" until the end of the tem¬ 
perature ramp. The same experiment is repeated four 
more times with fresh distilled water, and the results are 


pictured in Fig. [5jb)-(e) . The increase of G' is still vis¬ 
ible, but delayed from one experiment to the next and 
eventually no longer observed during the fifth experiment 
Pig - 1 e)]. Furthermore, comparing the surface state of 
the upper plate before the first experiment [Fig. [5][f )] and 
after the fifth experiment [Fig. [5j[g)], one observes an ox¬ 
ide layer visible by the naked eye in the latter case. The 
plate is indeed covered by some powdery material that 
can be removed easily with the finger which confirms the 
scenario discussed above: the corrosion of the duralumin 
plate nucleate aluminum oxide particles that increase the 
contact line elasticity leading to the premature increase 
of G'. Note that the bottom Peltier plate used in all 
the exper iments reported in the present article (except 
in section IIE1) is Teflon-coated which explains why the 
oxidation is only taking place at the rotor. 

In conclusion, one should be aware that sandblasting 
a duralumin plate makes it rough, but also promotes its 
slow oxidation by the agar sol, which leads to an artificial 
increase of G' that should not be interpreted as the early 
onset of the gel formation. We shall emphasize that such 
issue is absent when using commercial plates made of 
stainless steal. In section [TTTl all the gelation experiments 
are performed with an upper plate made of passivated 
duralumin to prevent any premature increase of G'. 


F. Solvent evaporation, oil coating and related artifacts 

Agar gels are mainly composed of water and evapora¬ 
tion has a strong impact on both the gelation dynamics 
and the gel linear properties. For gelation experiments 
performed under rapid cooling rates, i.e. T > l°C/min, 
evaporation is efficiently prevented by a solvent trap filled 
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FIG. 5. (Color online) (a)-(e) Temporal evolution of G' and G" of distilled water during a series of five consecutive experiments 
with mixed boundary conditions (rough duralumin upper plate and smooth bottom Peltier plate). Each experiment corresponds 
to a new water sample and consists in oscillations of small amplitude (7 = 1 %, / = 1 Hz) during a decreasing ramp of 
temperature from T = 70°C to 20°C at a cooling rate T = l°C/min with a solvent trap is filled with water. The zero normal 
force protocol is applied with F,y = (0.0 ±0.1) N and an initial gap value eo = 500 fim. Note that during each experiment, 
the gap remains constant equal to the initial value, which proves that there is no evaporation. The rotor has been sandblasted 
anew only before the first experiment displayed in (a). Pictures of the duralumin rotor: (f) before the first experiment and (g) 
after the fifth experiment. The scale is set by the rotor diameter of 40 mm. 


with water. However, for long lasting experiments con¬ 
ducted at lower cooling rates, evaporation becomes non 
negligible despite the solvent trap. A trick commonly re¬ 
ported in the literature consists in adding a thin layer of 
a non volatile oil at the periphery of the plate to isolate 
the sample from the ambiant ai j 20 * 22 * 6 ^ . The purpose 
of this section is to demonstrate that such method effi¬ 
ciently prevents the evaporation and to investigate quan¬ 
titatively the impact of the oil layer on the rheological 
measurements. 

To this end, we perform two gelation experiments of 
long duration at a cooling rate T = 0.2°C/min using 
a passivated sand-blasted duralumin rotor and the zero 
normal force protocol. The first gelation experiment is 
performed with water in the solvent trap [Fig. [6ja) and 
(b)], while in the second gelation experiment the sample 
is surrounded by a thin oil layer [Fig. |6j[c) and (d)]. In the 
first case, we observe visually that the water in the sol¬ 


vent trap topping the rotor has evaporated after about 
20 min, which coincides with the time above which G' 
shows non-negligible values [Fig. [6j[b)]. The elastic and 
viscous moduli intersect at t ~ 0.6 h while the gap starts 
decreasing. Both G' and G" slowly grow for about 2 h 
before showing a steep increase associated with the gela¬ 
tion. Non-negligible evaporation is further supported by 
two other observations. First, the large gap decrease of 
about 16% illustrated in Fig. [6](a) shows that the sample 
experiences water loss and that the increase of G' and G" 
between 0.5 and 2.5 h is artificial. Indeed, in the second 
experiment where the evaporation remains negligible be¬ 
cause of the oil surrounding the sample, the gap decreases 
by less than 1% over the entire experiment [Fig. [6j[d)] . 
Second, in the experiment conducted with water in the 
solvent trap, G' goes through a maximum at i ~ 4.2 h be¬ 
fore decreasing at larger times [Fig. [6j[b)— inset], whereas 
in the experiment using oil to prevent evaporation, G' ex- 
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FIG. 6. (Color online) Gelation dynamics of a 1.5% wt. agar solution induced by decreasing the temperature from T = 70°C 
to 20°C at a low cooling rate T = 0.2°C/min under a controlled normal force Fn = (0.0 ± 0.1) N. Temporal evolution of (a) 
the gap size e (initial value eg = 500 /im) and (b) the viscoelastic moduli G' and G" in semilogarithmic scale measured at 
/ = 1 Hz with the adapted strain protocol and an oscillatory frequency of / = 1 Hz. The inset is a zoom of G 1 vs time plotted 
in linear scale for 3.5 < t < 4.5 h. The solvent trap is filled with water, (c) and (d) same as in (a) and (b) except that the 
sample is surrounded by a thin oil layer. 


hibits a monotonic behavior and tends toward a constant 
value after 5 hours [Fig. |§d) -inset]. 

These results demonstrate quantitatively that for a low 
cooling rate of 0.2°C/min, the sample experiences ex¬ 
cessive evaporation when the solvent trap is filled with 
water. The use of an oil layer at the sample periph¬ 
ery prevents the solvent evaporation and allows one to 
access the steady state elastic properties of the geP-. 
However, if the use of oil leads to quantitative measure¬ 
ments in steady state, G' shows non negligible values 
from t ~ 0.85 h while the sample is still liquid, and long 
before the steep increase attributed to the gelation which 
takes place at f ~ 2.7 h [Fig. |6jd)] . Such behavior is 
reminiscent of the artifact described in the previous sub¬ 
section |lTE"2[ Nonetheless, the experiments reported in 
Fig-0 are conducted with a passivated duralumin rotor, 
thus the premature increase of G' is here related to an¬ 
other mechanism than the corrosion of the upper plate. 
We show in the remainder of this section that the prema¬ 
ture growth of G' is an artifact related to the presence of 
the surrounding oil layer. 

To study the effect of the surrounding oil alone, i.e. in 
the absence of the gelation, we perform a time sweep ex¬ 
periment on deionized water surrounded by a thin layer 
of oil. The experiment is conducted at constant temper¬ 
ature (T = 20° C) using a passivated duralumin upper 
plate and a transparent bottom plate made of PMMA 
(gap e = 500 /im). Although no gelation can take place, 
G' increases and becomes larger than G" from t > 1.6 h 


to reach a value G' — 10 Pa after six hours [Fig. [7j[a)]. 
Moreover, images taken through the transparent bottom 
plate shows that concomitantly to the growth of G', the 
oil rim slowly invades the gap at certain location along 
the sample periphery, and spreads on the surface of the 
upper metallic plate [Fig. [%J[b)-(f)]. Therefore, the pre¬ 
mature growth of G’ reported in Fig. id) during the 
gelation of the agar sol is most likley related to the inner 
migration of the oil/water interface. 

To further explore the role of the upper metallic 
plate on such an invasion process, we repeat the time 
sweep experiment on deionized water surrounded by a 
thin oil layer, with symmetric boundary conditions us¬ 
ing two identical transparent plates made of PMMA 
(constant temperature: T = 20°C and constant gap: 
e = 500 /mr). For more than 65 h, G' = 0 and only 
G" shows non-negligible values [Fig. |8j[a)] , while no mi¬ 
gration of the oil/water interface is observed by the naked 
eye [Fig.^b) -(e)]. These results confirm that the growth 
of G' and the migration of the oil/water interface ob¬ 
served in Fig. [7] are correlated. Furthermore it shows 
that the invasion of the gap by the outer oil layer, and 
therefore the premature increase of G', are directly re¬ 
lated to the use of metallic boundary conditions. 

Finally, if a more detailed analysis of the mechanism 
responsible for the oil invasion of the gap is out of the 
scope of the present publication, we discuss the weak 
linear increase of G" visible in Fig. |8j[a) which is instruc¬ 
tive. The increase of G" results from the slow pervapo- 
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FIG. 7. (Color online) (a) Evolution of the elastic and viscous moduli of deionized water surrounded by a thin layer of 
sunflower seed oil. Measurements are performed at T = 20°C in a plate-plate geometry with a constant gap e = 500 fi m, 
under small amplitude oscillatory shear ( 7 = 1 % and / = 1 Hz). The bottom plate is a transparent piece of PMMA, while 
the upper plate is made of sandblasted and passivated duralumin, (b)-(f) Snapshots of the periphery of the plate taken from 
below, through the transparent bottom plate at times t = 1, 2, 4, 6 and 8 h [these times are indicated by open symbols (□) in 
graph (a)]. The images show that the surrounding oil which is stained in red with a fat soluble dye (Sudan III), slowly invades 
the gap. 


ration of water through the oil layer. Indeed, a careful 
spatio-temporal analysis of the position of the oil/water 
interface pictured in Fig. [8])b)-(e) shows a net displace¬ 
ment SI of the oil/water interface after 65 hours® ranging 
between 30 and 50 /iin depending on the angular posi¬ 
tion along the interface [see Fig. 15 in the appendix]. 
Such a displacement corresponds to an average velocity 
of the oil/water interface of about 0.17 nm.s -1 . Assum¬ 
ing a near-diffusive transfer, the migration velocity v of 
the oil/water interface should scale as PD/l ;, where D 
is the diffusion coefficient of water in oil, P is the parti¬ 
tion coefficient of water at the oil/water interface and £ 
denotes the thickness of the oil layer®. For water perme¬ 
ation through a vegetable oil/water interfac^® one has 
P ~ 1.6 x 10 -3 and D ~ 10 -10 m 2 .s -17 h Considering 
an oil layer of thickness £ ~ 1 mm [Fig. [8jb)] , the above 
relation leads to v ~ 0.16 nm.s -1 in very good agreement 
with the experimental observations, which confirms that 
the pervaporation scenario accounts for the increase of 
G" in Fig.^a). 


In conclusion, using a layer of oil to surround the sam¬ 
ple efficiently prevents the evaporation and allows one 
to measure reliable steady-state values of the viscoelastic 
moduli. Nonetheless, one should keep in mind that the 
presence of an oil layer together with metallic bound¬ 
ary conditions may result in the slow migration of the 
oil/water interface at the periphery of the plates, in turn 
leading to the premature growth of G'. As a consequence, 
any early evolution of the elastic modulus during the for¬ 
mation of an agar(ose) gel monitored with a surrounding 
oil layer should be considered with caution. In particu¬ 
lar, the premature growth of G' reported in ref. 2 ^ in the 
early stage of a gelation experiment is likely to be artifi¬ 
cial and may not be the signature of the formation of an 


intermediate structure and/or a pre-gel inside the gap. 


III. IMPACT OF THERMAL HISTORY ON THE GEL 
MECHANICAL PROPERTIES 


In section[TTJ we have shown that the zero normal force 
protocol is an efficient way to monitor the heat-induced 
formation of agar gels. This method, which may be used 
with either rough or smooth boundary conditions, pre¬ 
vents the loss of contact between the gel and the plates 
and allows one to measure reliably the terminal values of 
the gel elastic and viscous moduli. In what follows, the 
zero normal force protocol is applied to determine the 
influence of the thermal history on the linear properties 
of a 1.5% wt. agar gel. We first quantify the impact of 
the cooling rate T on the gel elastic properties, before 
turning to the role of the final temperature T f reached 
at the end of the linear ramp, for a fixed cooling rate 
of l°C/min. These experiments are conducted with the 
mixed boundary conditions described in section IIB the 
upper duralumin plate is rough and passivated to prevent 
any oxidation (see subsection IIE 21, while the bottom 
Peltier plate is smooth and Teflon coated. Unless stated 
otherwise, the gelation experiments reported in this last 
section are conducted using an oil layer around the sam¬ 
ple to prevent evaporation (see section IIFI. 


A. Impact of the cooling rate 

In order to determine the impact of the cooling rate 
T on the gel elastic properties, we have performed a se¬ 
ries of gelation experiments on a 1.5% wt agar solution, 
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FIG. 8 . (Color online) (a) Temporal evolution of the elastic and viscous moduli of deionized water surrounded by a thin layer 
of sunflower seed oil. Measurements are performed at T = 20° C in a plate-plate geometry with a constant gap e = 500 pm, 
under small amplitude oscillatory shear (7 = 1 % and / = 1 Hz). The bottom and the upper plate are made of PMMA and 
transparent, (b)-(e) Snapshots of the periphery of the plate taken from below, through the transparent bottom plate at times 
t = 1, 24, 48 and 65 li [these times are indicated by open symbols (□) in graph (a)]. The surrounding oil is stained in red with 
a fat soluble dye (Sudan III). 


at different cooling rates ranging from 0.1 to 10°C/min. 
For each experiment, the temperature is decreased from 
70°C down to Tf = 20°C. The zero normal force proto¬ 
col is applied to monitor the gelation, and the value of 
the strain is further adapted to the value of the elastic 
modulus as discussed in IID to properly determine the 
intersection of G' and G". 


Results are displayed in Fig. [9] The formation dy¬ 
namics of the gel strongly depends on T: gelation oc¬ 
curs sooner for larger cooling rates. More quantitatively, 
the gelation time t gi defined here as the time at which 
G'{t g ) = G"(t g ) is inversely proportional to the cooling 


rate T and allows us to rescale the temporal evolution of 
the elastic modulus into a single mastercurve [Fig. |9][b)]. 
Such a rescaling hints at a unique gelation scenario that 
is independent of the cooling rate. The gelation temper¬ 
ature, defined as T g = T{t g ) decreases with the cooling 
rate: the gelation occurs at T = 37°C at T = 0.1°C/min, 
while the agar sol remains liquid down to 31° C for a cool¬ 
ing rate of 10°C/min Fig. |9jd)] . Interestingly, this evo¬ 
lution does not affect the steady-state properties of the 
gel at T = 20° C since the terminal values G'j and G'j are 


independent of T, within error bars [Fig. |9|[c)]. Scanning 
cryoelectron microscopy (cryo-SEM) shows that indeed 
the microstructure of gels prepared at different cooling 
rates is qualitatively the same (Fig. 101 which supports 
the idea that lower cooling rates only delay the gelation, 
but do not affect the terminal elastic properties of a 1.5% 
wt. agar gel. 


This observation strikingly contrasts with results from 
the literature where agarose gels prepared at con- 
stan t gap width appear stronger for decreasing cooling 
rates^ 7 H We attribute such a discrepancy to the par¬ 
tial loss of contact between the gel and the geometry in 


previous studies, as already discussed in Fig. [lja) and 
(d). Indeed, the formation of agar gels causes the sample 
contraction which is robustly observed to become more 
important for increasing cooling rates, independently of 
the presence of an oil rim [Fig. [9jd)]. The sample con¬ 
traction increases from about 1% for T < l°C/min, up 
to 5%, at T = 10°C/min. In the latter case, the rapid 
cooling freezes the microstructure far from its minimum 
energy state and the relaxation of internal elastic stresses 
then favors the contraction of the network. For a fixed 
gap, the contraction is sufficient to trigger a partial loss 
of contact between the gel and the geometry leading to 
a lower estimate of the elastic modulus, exactly as evi¬ 
denced in the experiment reported in Fig. [ija) and (d). 
Besides, the present observation showing that larger cool¬ 
ing rates lead to larger contractions is compatible with 
the apparent decrease of the elastic modulus reported 
for increasing cooling rates in experiments previously re¬ 
ported in the literature at constant gap width. Hence 
the relevance of the zero normal force protocol, espe¬ 
cially for large cooling rates where the gel contraction 
becomes non-negligible. Finally, one should keep in mind 
that the value of T selects the gelation temperature T g 
(roughly estimated from the intersection of G' and G"), 
and does not affect the gel linear rheology determined 
at 20° C, at least over two decades of cooling rates from 
0.1 to 10°C/min. These two results are not incompat¬ 
ible. Indeed, the elastic modulus of the agar gel scales 
as a power-law of the concentration in agarose, with an 
exponent close to 2.2 25 . Therefore a 5% increase in the 
agar concentration due to the gap decrease during a rapid 
cooling results in a 12% increase in G'. The uncertainty 
in the measurement of G' is of the same order of mag¬ 
nitude as emphasized in section m which explains why 
G' is reported as roughly constant without significant de- 
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FIG. 9. (Color online) (a) Evolution of the elastic modulus G' vs. time t when decreasing the temperature of a 1.5% wt. 
agar solution from 70°C down to 20°C. Colors from black to yellow encode different cooling rates: T = 10, 5, 2, 1, 0.5, 0.2, 
0.1°C/min. (b) Same data set plotted vs. the normalized time t/t g , where t g is the gelation time defined by the intersection 
of G' and G". Inset: gelation time t g vs. the cooling rate T. The black line corresponds to the best linear fit of the data in 
logarithmic scale, which equation is: t g T = 34.9°C. (c) Terminal values G'f (■) and G'j (A) of the elastic and viscous moduli 
respectively vs. the cooling rate T. The horizontal dashed lines stand for the mean values of the gel viscoelastic moduli: 
G'f = (24.2 ±0.6) kPa and G'j = (340 ± 10) Pa. For the sake of clarity, the typical error bar which corresponds to the dispersion 
of the results obtained by repeating the experiment a few times is indicated only on one data point for G'f. (d) Gelation 
temperature T g vs. the cooling rate T. (e) Relative variation of the gap thickness e vs. the cooling rate T. Each experiment 
is performed with a thin layer of oil surrounding the sample. Open squares (□) denote experiments performed with water in 
the solvent trap and no surrounding oil layer. 


pendence upon the cooling rate. 


B. Impact of the temperature drop 

The cooling rate is now fixed to T = l°C/min and ex¬ 
periments are performed by varying the final temperature 
Tf at the end of the ramp so as to quantify the impact of 
the temperature drop on the gel mechanical properties. 
The sample is loaded as a liquid into the preheated gap of 
the rheometer, and the temperature is swept from 70°C 
down to Tf, where Tf is chosen between 20 and 38°C. 
As illustrated in Fig. ITlTa) , the terminal temperature Tf 
impacts both the gelation dynamics and the steady-state 
value of the elastic modulus. 

Imposing a cooling rate T = l°C/min, fixes the gela¬ 
tion temperature to T g ~ 35.5°C [Fig. [9](d)]. We will 
therefore discuss the gelation dynamics whether Tf is 
larger or smaller than 35.5°C. For Tf < 35.5°C, the gela¬ 
tion always starts at T g = 35.5°C confirming that the 
intersection of G' and G" mainly depends on the cooling 


rate [Fig. [lTJb)]. Moreover, the time tf for the elas¬ 
tic modulus to reach a steady-state value (within 5%) 
increases with Tf, which accounts for the increasing du¬ 
ration of the gelation as the terminal temperature is in¬ 
creased [Fig. [Tljc)]. Surprisingly, for Tf > 35.5°C, the 
gelation dynamics becomes slower but the gelation is still 
observed, at least up to T = 38°C. Within this temper¬ 
ature range, the slowdown of the dynamics now results 
from the increase of both the gelation time t g and the 
time tf to reach the steady state, i.e. for increasing val¬ 
ues of T f , the intersection of G ' and G" occurs later 
and once G' > G" it takes more time for G 1 to reach 
G'f [Fig. 11 b) and (c)]. Finally, for increasing terminal 


temperatures Tf one observes the formation of gels with 
decreasing elastic modulus G'f [Fig. fTTVd)] . Such a be¬ 
havior is observed whether Tf is smaller or larger than 
T s (l°C/min) ~ 35.5°C, which suggests that the gel for¬ 
mation scenario is likely to be the same on both sides of 
T 0 (l°C/min) although the dynamics is strongly slowed 
down for T f > T fl (l°C/min). 


Indeed, here again the various gelation curves G'(t ) 
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FIG. 10. Cryo-SEM images of 1.5% wt. agar gels prepared at different cooling rates: (a) T = 0.2°C/min, (b) T = l°C/min, 
and (c) T = 10°C/min. In each case, a gel sample is mounted on a pin stub and immersed in liquid nitrogen for about 5 min, 
before being placed in the low vacuum preparation chamber of a scanning electron microscope (JEOL 6700F). The chamber is 
cooled down to T = — 90°C and the sample cut in situ with a scalpel. The temperature is increased at a rate of 5°C/min up to 
T = — 50° C and maintained constant for about 5 min to sublime the water frozen inside the gel. The temperature is decreased 
back to T = —85°C. The sample is coated with a nanolayer of gold-palladium, before being cooled down to T = — 160°C and 
finaly introduced in the observation chamber of the microscope. Samples are imaged in SEI mode, at 5 kV. 


obtained for different terminal temperature Tf can be 
rescaled into a single master curve by plotting G'/G'f 
versus (t — t g )/(tf — t g ) [Fig. [TT)(e)]. The latter function 
provides the best rescaling among other temporal func¬ 
tions such as t/tf, t/t g , etc. which demonstrates that 
the duration tf — t g separating the beginning of the gela¬ 
tion from the steady state is the relevant unit of time 
associated with the gelation of 1.5% wt. agar gels. Fur¬ 
thermore, the existence of a master curve confirms the 
existence of a single gelation scenario, that is suppos¬ 
edly the same as the one evidenced by the experiments 
performed at different cooling rates (section III AI. This 
result also indicates that the gel elastic modulus is mostly 
a function of the terminal temperature Tf reached at the 
end of the temperature ramp, and that G'f(Tf) reflects 
the temperature dependence of a gel with a single mi¬ 
crostructure, whether it has been formed above or below 
35.5°C. This point is discussed in more detail in the fol¬ 
lowing paragraph. 


C. Impact of a temperature plateau during cooling 

To determine whether the gels that are obtained 
at different terminal temperatures on both sides of 
T s (l°C/min) ~ 35.5°C display fundamentally different 
properties, we perform a last series of experiments. Gela¬ 
tion experiments of 1.5% wt. agar solutions similar to the 
ones reported in section |HI B| are performed by decreasing 
the temperature at l°C/min from 70°C down to T = T p , 
where T p is chosen between 20 and 38° C and maintained 
for 8 h before resuming the cooling down to Tf = 20° C. 
The result of a temperature stop at T p = 33° C is given in 
Fig^a). The elastic modulus reaches a constant value 
while T = T p in agreement with the data reported in 
Fig. m The subsequent drop of temperature down to Tf 
leads to a novel increase of G' which terminal value, G/, 


is compatible with the value measured by direct cooling 
from 70°C to 20°C, i.e. without any plateau of temper¬ 
ature. The experiments repeated at different values of 
T p further support this conclusion: the terminal values 
reached by G' and G" are independent of T p [Fig. 121(b)], 
whether the plateau of temperature is performed above 
or below the gelation temperature Tf = 35.5°C selected 
by the cooling rate. Such values of G'f and G'/ are also 
compatible within error bars with the ones reported in 
Fig. [|(c) for different cooling rates. This result shows 
that the gelation of a 1.5% wt. agar solution is mainly 
controlled by the spinodal demixing expected to take 
place over shorter timescales than the timescale associ¬ 
ated with the cooling process (i.e. 50 min). As a result, 
the elastic and viscous moduli of the gel formed upon 
cooling are fixed by the terminal temperature Tf and 
the thermal history poorly impacts the gel linear prop¬ 
erties. Finally, the gel contraction slightly varies from 
one experiment to the next, but does not show any sig¬ 
nificant trend with the temperature of the plateau, rang¬ 
ing between 0 and 2%, with a mean value Ae/e ~ 1 % 
[Fig.[l2|(c)]. Again, such values justify the use of the zero 
normal force protocol instead of a constant gap width to 
determine the gel elastic properties. 


IV. DISCUSSION AND CONCLUSIONS 

In the first part of this work we have shown quantita¬ 
tively that the zero normal force protocol is more adapted 
to monitor the gelation dynamics of agar solutions than 
the traditional constant gap width procedure. Indeed, 
the gelation of a 1.5% wt. agar solution results in the 
sample contraction as evidenced by the large negative 
normal forces measured in experiments performed at con¬ 
stant gap width. Such a contraction in a cell of constant 
volume leads to the strain hardening of the gel, which 
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FIG. 11. (Color online) (a) Evolution of the elastic modulus G' vs. time t while cooling the agar solution at a constant rate 
of T = l°C/min, from T = 70°C down to Tf. The colors, from yellow to black, stand for different end-of-ramp temperatures: 
Tf = 20, 30, 32, 34, 35, 36, 37 and 38°C. (b) Gelation temperature T g as determined by the intersection of G' and G " vs. the 
terminal temperature Tf. (c) Time tf for the elastic modulus to reach the steady state value vs. the terminal temperature 
Tf in semilogarithmic scale, (d) Steady-state value of the elastic modulus G'f vs. Tf. (e) Normalized elastic modulus G'/G'f, 
where G'f stands for the steady state value of G 1 , vs. (t — t g )/(tf — t g ). Each experiment is performed with a thin layer of oil 
surrounding the sample. 


explains the artificial drift of the gel elastic modulus to¬ 
wards ever larger values that is commonly encountered in 
the literature, even after the gelation is over. Although 
supplemental observations of the gel formation in a ge¬ 
ometry of constant gap width are needed to fully eluci¬ 
date the local scenario associated with strain hardening 
and the gel debonding from the walls [see for instance 
re fP), we have shown that the zero normal force proto¬ 
col allows one to determine more reliably the steady-state 
values of the gel elastic properties. Furthermore, we have 
demonstrated that the zero normal force protocol makes 
it possible to monitor the gelation with smooth boundary 
conditions. Therefore, such protocol ensures a good ther¬ 
mal contact between the sample and the bare metallic 
geometry compared to constant gap width experiments 
which require the use of sandpaper (or an equivalent coat¬ 
ing) to prevent wall slip and/or the gel debonding from 
the walls. As such, the zero normal force protocol ap¬ 
pears very promising to investigate more quantitatively 
the gelation or the solidification dynamics of samples dis¬ 
playing large volume changes during the sol-gel or the 
liquid-solid transition. In particular, we anticipate that 
the zero normal force protocol will find relevant applica¬ 
tions for cocoa butter rheology' 4 76 and mo re gen erally 
in the rapidly growing field of hydrate slurries'' 78 where 
rheological measurements appear extremely delicate due 
potentially to the relative dilation of the sample com¬ 


pared to that of the geometry under complex thermal 
protocols. 

Moreover, we have unraveled two subtle artifacts re¬ 
lated to the contact of an aqueous solution with metal¬ 
lic boundary conditions that may lead to a premature 
increase of the elastic modulus long before the start of 
the gelation. On the one hand, the corrosion of dura¬ 
lumin plates by the oxidative agar sol may enrich the 
air/solution or the oil/solution interface with nanopar¬ 
ticles. Such nanoparticles increase the elasticity of the 
contact line before the start of the gelation, leading to an 
artificial increase of the elastic modulus, with G' > G". 
Note that the latter issue is specific to duralumin and is 
not observed for plates made of stainless steel. On the 
other hand, the oil layer traditionally added around the 
sample to prevent evaporation may also affect the rheo¬ 
logical measurements while the sample is still liquid. The 
oil may slowly invade the gap leading to a premature in¬ 
crease of G'. The role of the plates wetting properties on 
the oil invasion scenario remains to be investigated. 

Finally, in the second part of this article, we have ap¬ 
plied the zero normal force protocol to determine the in¬ 
fluence of thermal history on the linear rheology of a 1.5% 
wt. agar gel containing 1% wt. of agarose. Our study 
demonstrates that the gel elastic properties are mostly 
controlled by the terminal temperature at the end of the 
cooling process, and that contrarily to gels with higher 
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FIG. 12. (Color online) (a) Evolution of the elastic modulus 
G' (solid line, left axis) and of the temperature T (dashed 
line, right axis) versus time t. The temperature is decreased 
at a cooling rate T = l°C/min, from T = 70°C down to 
T = T p = 33° C. The latter temperature is maintained for 
8 h before resuming the cooling down to Tj = 20° C at 
T = l°C/min. (b) Terminal values G'f (■) and G'f (a) of 
the elastic and viscous moduli vs the temperature T p of the 
plateau. For the sake of clarity, the typical error bar deter¬ 
mined by repeating the experiment is indicated only on one 
data point. The horizontal dashed lines stand for the mean 
values of the gel viscoelastic moduli: G'f = (23.2 ± 0.5) kPa 
and G'f = (350 ± 10) Pa. (c) Relative variation of the gap 
vs. T p . Each experiment is performed with a thin layer of oil 
surrounding the sample. 


agarose contenP^, neither the cooling rate nor the ther¬ 
mal path affect the gel mechanical properties. Such a 
result suggests that the local gelation scenario is mostly 
controlled by the rapid spinodal demixing of the agar so¬ 
lution in contrast with the complex formation scenario 
reported irPH on a 2% wt agarose gel. As a practical con¬ 
sequence, the thermal parameters allow one to modulate 


the gel formation dynamics. Namely, the crossover be¬ 
tween the elastic and viscous modulus is selected by the 
cooling rate T, while for a given value of T, the choice 
of the terminal temperature relative to the temperature 
at which the crossover of G' and G" takes place makes it 
possible to strongly slow down the gelation. These results 
should help in the optimization of manufacturing lines of 
agar-based products such as Petri dishes or specific elec¬ 
trophoretic gels, for which the cooling rate and the delay 
to complete gelation are the most constraining factors 
regarding the production rate. More generally, the zero 
normal force procedure should help performing more reli¬ 
able rheological investigation of non-isochoric processes. 
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FIG. 13. (Color online) (a) Temporal response of the gel elas¬ 
tic modulus G' to an increase of the gap width of Ae = 6 /rm 
triggered at t = 0.5 h and suppressed after 1 h. Ae denotes 
the excess gap width above the minimum value reached at 
t=0.5h, at the end of the gelation process under controlled 
normal force, and taken as an arbitrary reference, (b) Tem¬ 
poral evolution of the normal force Fn concomitant to the gel 
response. 
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FIG. 14. (Color online) Temporal evolution of the normal force Fjv and the gap width e (a)-(c) together with the elastic and 
viscous moduli (d)-(e) during the gelation of a 1.5% wt. agar solution induced by decreasing the temperature from T = 70° C 
to 20°C at a cooling rate T = l°C/min. The gelation experiment is repeated three times using the zero normal force protocol 
and applying different constant strain amplitudes: 7 = 0 . 01 %, 0 . 1 % and 1 % in respectively the first, second and third column. 
Terminal values of the gel elastic modulus are respectively: 22.3, 21.8 and 22.3 kPa, and compatible within error bars. For all 
the experiments the oscillation frequency is 1Hz. (g)-(i) Strain sweeps experiments performed at 1Hz on the three gels obtained 
in (d), (e) and (f) after two hours, once they are fully gelihed. The vertical dashed line denotes the critical strain 7 C above 
which G' decreases. 


APPENDIX 

The first part of the appendix provides a proof that the 
slow increase of the elastic modulus reported in Fig. [lje) 
during a gelation experiment performed at constant gap 
width, is indeed due to the strain hardening of the sam¬ 
ple. We have performed the following supplemental test: 
an agar gel is prepared between parallel plates (initial 
gap eo =500 grn) using the zero normal force protocol, 
by decreasing the temperature from 70°C down to 20°C 
at l°C/min. At the end of the gelation process, the gap 
has decreased by 1 % and the normal force is still equal to 
zero, in agreement with the results reported in Fig. [ljc) 
and (f) in the main text. We then switch the rheometer 
to a controlled gap mode and increase the gap width by 
Ae = 6 /iin at 0.1 /irn/s, which corresponds to a verti¬ 
cal extension of the gel of about 1%. In response, the 
gel elastic modulus reported in Fig. 13 ’a) shows a brutal 


increase of about 10% from 23 kPa to 25 kPa, while the 
normal force concomitantly becomes negative and relaxes 
towards F N ~ — 3N [Fig. 


13 


a)] which indicates that the 
gel is pulling on the upper plate. These results show that 


an increase of the gap width triggers an increase of the gel 
elastic modulus: the gel experiences strain hardening un¬ 
der extension. We thus conclude that the contraction of 
the sample at constant gap width also leads to the strain 
hardening of the sample, which explains the ever increas¬ 
ing values of G' reported during gelations performed at 
constant gap width [Fig. [l|e)]. Moreover, when the gap 
is brought back to its initial position after 1 h, the normal 
force relaxes back to zero, while G' decreases to a value 
slightly lower than the initial one (Fig. [l3]) . Such a dis¬ 
crepancy is likely due to a partial loss of contact between 
the gel and the plates, but remains to be investigated in 
more details. 

The second part of the appendix is dedicated to the 
impact of the strain amplitude during the gelation ex¬ 
periments. We have repeated the gelation experiment 
reported in the main text [Fig. 0 f)] at different strain 
amplitudes of constant value: 7 = 0 . 01 , 0.1 and 1 % while 
applying the zero normal force protocol. The results dis¬ 
played in Fig. |14| show quantitatively the same terminal 
values of the elastic modulus [G^ ~ 22 kPa in Fig. 

(f)] and the same sample contraction [Ae/e 
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~ 1 % 
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FIG. 15. (Color online) (a) Picture of the parallel plate geometry shot from below. Both the upper and bottom plates are 
made of transparent PMMA. The gap is filled with distilled water and surrounded by an oil rim of about 1 mm thick and 
stained in red with a fat soluble dye. The scale is set by the plate diameter of 40 mm. (b) Spatiotemporal diagram of the 
oil/water interface as a function of time and computed in the region of interest (ROI) emphasized by the yellow rectangle in (a), 
(c) Radial velocity v of the water/oil interface computed over 65 h vs. the angular position 9 of the ROI. The ROI reference 
angular position (9 = 0°) is indicated in (a) and the angle 9 is counted positively in the counterclockwise direction. 


Fig. [HJa) -(c)]. Therefore, as long as the strain ampli¬ 
tude remains lower than 1 %, the stress applied during 
the gelation performed with the zero normal force proto¬ 
col does not impact the gel formation dynamics nor the 
steady-state parameters of the gel. Yet the strain applied 
during the gelation does matter. Indeed, we have per¬ 
formed a strain sweep on the three gels which have been 
exposed to different oscillatory amplitudes during their 
gelation [Fig. [I4|(g)-(i)]. For all the experiments, we ob¬ 
serve that G’ decreases above a critical strain y c which is 
interpreted as the progressive debonding of the gel from 
the plate, as the gel remains intact and no macroscopic 
crack are visible at the end of the strain sweep. First, 
the critical strain y c above which G' decreases as a re¬ 
sult of the gel debonding is always larger than 1%. This 
result further justifies the strain adapted protocol intro¬ 
duced in section m of the main text, which consists 
in progressively decreasing the strain amplitude applied 
to measure the viscoelastic modulus during the gelation. 
Second and as a remarkable result, the critical strain y c 
above which G' decreases depends on the strain ampli¬ 
tude applied during the gelation. Gels which formation 
has been monitored with a smaller amplitude, maintain 
contact with the plates up to a larger value of 7 C . This 
observation indicates that small strain amplitude as low 
as 0 . 01 % during the gelation are required for an optimal 
adhesion between the sample and the plates. Finally, the 
evolution of the viscoelastic modulus in the vicinity of 
the critical deformation y c strongly differs depending on 


the strain amplitude applied during the gel formation. 
The gel which has been prepared at 7 = 0.01% shows 
a brutal drop in G' during the strain sweep [Fig. 14 (g)] 
which hints at an abrupt detachment from the plates. 
Gels which formation have been monitored under larger 
strains (7 = 0.1 and 1 %) display a much smoother de¬ 
crease in G' above y c [Fig. [±I h ) and (i)] which could 
correspond to a more progressive contact loss from the 
plate. This last point certainly deserves more investiga¬ 
tion. 


The third and last part of the appendix concerns the 
displacement of the peripheral oil/water interface sand¬ 
wiched between two transparent PMMA plates, as illus¬ 
trated in Fig. [ 8 ] in the main text. Let us recall that the 
gap is filled with distilled water and that the surrounding 
oil rim is stained in red with a fat soluble dye. We use a 
spatiotemporal analysis to measure the velocity of the in¬ 
terface. A picture of the entire geometry shot from below 
is reported in Fig. [l5ja). A rectangular region of interest 
(1.45 mmxl.2 mm, i.e. 39x32 pixels 2 ) is pictured as a 
yellow rectangle in Fig. [H^a) and defines the reference 
angular position 0 = 0 ° with respect to the plate rotation 
axis. The temporal evolution of the gray level along the 
oil/water interface is plotted as a spatiotemporal diagram 
in Fig. |15[ b) which allows us to measure the inclination 
of the grey level texture through a sub-pixel digital image 
correlation and detect a displacement of about one pixel 
over a duration of 65 h! The latter analysis applied every 
10 ° along the periphery of the oil/water interface reveals 
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that the oil moves systematically towards the rotation 
axis at a velocity v which depends on the angular posi¬ 
tion 0 and ranges between 0.4 and 0.8 /rni/h [Fig. |T5^c)] . 
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